Determination of Boundary Scattering, Intermagnon Scattering, and 
the Haldane Gap in Heisenberg Chains 
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Low-lying magnon dispersion in a S = 1 Heisenberg antiferromagnetic (AF) chain is analyzed 
using the non-Abelian DMRG method. The scattering length Ob of the boundary coupling and 
the inter-magnon scattering length a are determined. The scattering length Ob is found to exhibit 
a characteristic diverging behavior at the crossover point. In contrast, the Haldane gap A, the 
magnon velocity v, and a remain constant at the crossover. Our method allowed estimation of the 
gap of the S — 2 AF chain to be A = 0.0891623(9) using a chain length longer than the correlation 
length ^. 
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I. INTRODUCTION 

To form a better understanding of interacting many- 
body systems, it is very important to determine an effec- 
tive field theory and to clarify the low-energy physics 
involved. In the physics of low-dimensional quan- 
tum systems, considerable attention has been paid to 
the one-dimensional antiferromagnetic (AF) integer-spin 
, Heisenberg model following the discovery of the Haldane 
^ I gapj ^^'^^ Precise determination of the gap has been re- 
ported by several authors. Its massive elementary 
excitation, i.e., the magnon, has a relativistic dispersion 
relation, which is often described by a non-linear sigma 
model (NLSM),M>25i4i^ 

In particular, the S* = 1 AF Heisenberg chain has 
been widely studied both theoretically and experimen- 
tally. When open boundary conditions (OBC) are ap- 
plied to a S" = 1 AF chain, owing to the unique effective 
S = 1/2 spins at the ends, quasi-degeneracy appears be- 
tween the singlet ground state and a low-lying triplet 
state4^ Various attempts at boundary tuning, as ex- 
emplified by attachment of real S = 1/2 spins to maintain 
the high accuracy of the density matrix renormalization 
group (DMRG) method^^i^ have shown that deforma- 
tion of the boundary conditions can selectively modify 
the magnon wavefunction while maintaining the unifor- 
mity of the ground state. 

To form a better understanding of the physics involved 
in a finite chain under OBC, we can use the NLSM 
to describe the low-lying energy dispersion. Lou et al. 
have proposed usage of a form, \/A2 + sin^ fccff , for 
low-lying magnon dispersions, where fccff is the effective 
wavenumber.^^ Here, A denotes the Haldane gap and 
V is the velocity of the quasi particle. They described 
the asymptotic effects of boundary scattering and inter- 
magnon interactions in terms of the scattering lengths, 
Ob and a, which appear in fceff. When boundary tuning 
is applied by introducing an antiferromagnetic coupling 
Jcnd between the 5 = 1 spin chain and the extra real 
5 = 1/2 spin, these scattering lengths might be effected. 
This idea motivated us to study low-lying elementary ex- 



citations using both the DMRG and NLSM methods by 
describing the bulk properties and the boundary scatter- 
ing effects in terms of an effective theory. In this work, 
using the DMRG method, the energy dispersion of var- 
ious magnon modes was determined for 5 = 1 Heisen- 
berg systems with up to 2048 spins. Finite-size scaling 
analysis was performed to determine the boundary scat- 
tering length and the inter-magnon scattering length, in 
addition to A and v in the thermodynamic limit. We 
used a relation of the correlation length ^ ~ w/A, which 
is known to hold approximately in the integer-spin AF 
Heisenberg chain. ''^•^^ We found that changed sign 
around a critical value of Jcnd- This value should be 
identical to that required to make local quantities such 
as the local bond energy of the ground state and the spin 
density of long- wavelength magnons uniformi ^^i'^^i^°i^^ In 
addition, a divergence- like behavior of Ob was detected 
around this critical value denoted as J^^^. However, the 
inter-magnon scattering length was found to be constant 
at a = —0.383(6)^ irrespective of Jcnd- In this deriva- 
tion, A, V, and ^ were confirmed to be always indepen- 
dent of Jcnd in the thermodynamic limit. This allows 
the low-lying elementary excitations to be effectively de- 
scribed. The results indicated the presence of both itin- 
erating magnons (IMs) and boundary magnons (BMs) 
bound at the ends. At J^.^^, the diagonal magnetization 
induced by an IM shows a flat structure around the cen- 
ter of the system when L ^ ^, with L being the number 
of 5 = 1 spins. Both the diverging behavior of ab and 
the uniform distribution of the long- wavelength magnons 
confirm the realization of bulk characteristics in an ele- 
mentary excitation at the critical point J^^^, where the 
ground state also has a uniform nature around the center 
of the system. 

Furthermore, this work clearly resolves the problem 
pointed out by Todo and KatoK; there is disagreement 
between the DMRG'^^ and quantum Monte Carlo (QMC) 
simulation results'^^ with respect to estimation of the ex- 
citation gap in the 5 = 2 AF Heisenberg model. The 
reason for this disagreement might be an inappropriate 
scaling assumption in the DMRG study. This work ap- 
plies finite-size scaling analysis to the excitation gap in 
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the S = 2 AF chain, and shows for the first time that the 
corrected gap is within the error bar of the QMC value. 



II. EFFECTIVE HAMILTONIAN 

We consider a. S = 1 AF chain with boundary S = 1/2 
spins Sj with j = or L + 1 , which is described by the 
foUowing Hamiltonian. 
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S; • Si+i + Jcnd(so • Si + Sl ■ Sl+i), (1) 



where Si represents the S = 1 operator at the z-th site. 

The low-energy physics of the Hamiltonian in Eq. ([T]) 
can be understood by using an approximate mapping 
onto the NLSMi^i^iiiii^ We let L ^ oo, keeping X = Lb 
constant, with b being the lattice spacing. Taking into 
account the effective S' = 1/2 boundary modes s^^, we 
obtain the following expression, 

H,s = i?NLSM + A, [0(0) • sf + (-1)^0(A) • sf ] 
+A„ [1(0) ■ sf + 1(A) . sf ] 
+A', [ckiO) ■ So + (-1)^0(A) • SL+i] 
+A: [1(0) • So + 1(A) • SL+i] 
+ Jifd [sf • so + sf . Si+i] , (2) 

with the bulk part of the NLSM expressed as 



^ Jo 



1 

g \dx 



(3) 



where <p and 1 = {l/vg)(p x dt4> are low-energy Fourier 
modes of the spin operators with wave vectors near tt and 
0. The coupling parameter and the velocity are given 
= 25*. Since all the bare couplings are 
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antiferromagnetic, solutions for the bulk fields follow the 
Neumann boundary conditions (NBC) : d(f>/dx\x=o.x ~ 
0^ The A„ and A^ terms produce an effective boundary 
repulsive potential on an IM, and J^rfd 

is a renormalized 

coupling constant. 

The validity of this description is also confirmed by 
examining the spin density of an IM shown in Fig. [T] 
When Jcnd is larger than JcndJ the lowest triplet mode 
has itinerating behavior. Indeed, we see that (S^) ex- 
hibits a cosine- like behavior for Jcnd = 1.0 owing to both 
strong repulsive coupling via A„ and A^, and the NBC 
on (f){x). When J^nd approaches J^^^ ~ 0.51, the IM 
mode becomes uniform around the center of the chain 
but (Sf) exhibits damped oscillations near the two ends. 
This known solution suggests that the mode should con- 
tinuously change into an end mode sf in the low-energy 
eigenstate when Joff < Jcs- 

Thus, the dispersion relation for N itinerating magnon 
modes at low energies in the dilute limit may be simply 
reproduced by a nonrelativistic effective Hamiltonian for 
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FIG. 1. (Color online) Distribution of local magnetization 
(Sj) for a single magnon state with Stot ~ 1 for various values 

of Jcnd. 



N virtual particles, 
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HMnd) = ^ H ^ + H Vi^^ - ^j) 
AT 

+ 5I[^b(Jc„d,a;0 Vfc(Jcnd, A - Xi% (4) 
1=1 

where < < A, with a wavefunction obeying the Neu- 
mann boundary condition: djip{xi, ■ ■ ■ ,a;jv)Uj=o.x — 0. 
Here we use the Einstein relation m = A/v^. Effective 
short-range interactions between IMs and between an IM 
and a BM are represented by V and Vb, respectively. We 
expect they are short-range function with range of the 
order of the correlation length f . All of the effects of 
Jcnd are produced by the boundary potential Vb(Jcnd7 2;) 
In the asymptotic region, the effects of V and Vb appear 
as scattering phase shifts, which are represented by a and 

We now identify low-lying magnon modes. Each mode 
is specified by a total spin of S'tot- When Jcnd is small 
and positive, since we have two effective = 1/2 spins 
creating the bulk low- lying triplet and two real = 1/2 
spins, we need to polarize these four spins before we can 
create one IM. In this case, the effective chain length for 
the IM becomes L — 2ab and k^s — 7r/(_L — 2ab), when 
the system is about two times longer than the correlation 
length ^. Therefore, we have the relation: 



£^32 — 



sin^ 



L — 2ab 



(5) 



where Eji = Ej — Ei and Ej and Ei are the lowest energy 
of the S'tot — j and S'tot — i states. The energy spectrum 
i?42 for two IMs is given by 
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L — 2ab 



(6) 
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where we use the small-fc approximation for the magnon- 
magnon phase shift. When Jcnd becomes large enough, 
the effective boundary 5 = 1/2 modes couple strongly 
with the real S — 1/2 spins. In this condition, the low- 
lying magnon states are IMs, and the formulas for ii'io 
and £^20 are, respectively, similar to Eq. ^ and Eq. 
©. Thus, we can conclude that a crossover value of J^^^ 
exists, where the low energy spectrum changes qualita- 
tively. 



III. NUMERICAL RESULTS 

We used the non-Abehan DMRG method (NA- 
DMRG)^^ to estimate the energy spectrum of the low- 
est ^tot = 0,1,2,3, and 4 states for finite systems. 
Numerical convergence during finite system sweeping 
was accelerated by the use of a wave function predic- 
tion method. ^^"^^ Since the number of kept states for 
the block spin is up to = 512, the truncation er- 
ror is smaller than 1.0 x 10^^ in the lowest 5tot = 4 
state. This corresponds to a number of kept states of 
■nis^ ~ 2500 - 2700 in the standard DMRG. In this case, 
the numerical cost of the standard DMRG is about 110 
- 140 times higher than that of NA-DMRG, because in 
the DMRG it varies as the cube of the number of kept 
states. The system size L -|- 2 is up to 2048, where the 
two extra spins indicate the boundary 5=1/2 spins. 

The energy of a single IM as a function of the system 
size is shown in Fig. [51 The target energy spectrum is 
£■32 when Jend = 0, and Eiq when Jend = 0.6 or 1. To 
estimate A, v and Cb, we generated sequences yl*(Lo + 2) 
for different values oi L = Lq, where A*{L) denotes finite 
values oi A — A,v and Ob in the thermodynamic limit. 
The sequences were determined by least square fitting 

]^^2at, energies 



with the function a/A'^ + v'^ sin^ 



of L 2 = 2*^(Lo + 2), where ^ = 0, ±1. The value of A 
was estimated by power-law extrapolation with elements 
of A* (512) and A* (1024). The estimation error was taken 
to be \A — A* (1024) I . Based on the optimum boundary 
scattering length ai,{ Jcnd) for each Jend, we found a uni- 
versal finite size dependence for a fixed energy gap A 
and spin velocity v. As a result, we showed that only the 
boundary scattering length ab was affected by changing 
Jcnd J whereas A and v were independent of Jcnd (See Ta- 
ble |l|. This result is consistent with the effective model 
in Eq. (g]). 

The estimated values of A, v, and ^ = v/A are con- 
sistent to within A = 0.4104792485(4), v = 2.46685(2) 
and ^ = 6.00967(5), respectively, except for a some- 
what larger error at Jcnd = 0.4 and 0.6, which are 
closest to J^nd- Since our data is obtained by extrap- 
olation using system sizes larger than those treated in 
former studies,—!^ our results show meaningful differ- 
ences. The reported value of ab(Jend = 0) = — 1 in 
Ref. f^] is about three times larger than our result of 
ab(Jend = 0) = -0.3748(1). The value of ab(Jond) 
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FIG. 2. (Color online) Single magnon energy with Jcnd un- 
der the condition = 512 (m^z ~ 2500). The dotted line 
represents sin^(7r/(L — 2ab)). 



changes rather dramatically with Jcnd, (Fig- El with a 
change in sign even occurring around Jcnd ^ 0.5. The 
values seem to diverge around Jpnd- When the boundary 
scattering length becomes ab —00, fccfr = Tr/{L — 2ab) 
approaches zero, and the energy of the lowest IM is al- 
most at its minimum value. A, and is independent of 
Jcnd- This is consistent with the report in Ref. pfi'Soj^ 
However, we should note that the above picture holds 
only when L 3> flb, requiring a high performance simu- 
lation tool such as NA-DMRG. 

In the same manner, using the estimated A, v, and 
flb, we determined the inter-magnon scattering length 
a. The target energy spectrum is E42 when Jcnd < 
J^jj(j, and is E20 when Jend > ^^end- With a com- 
mon a, universal behavior is observed in the large L re- 
gion. The estimated values of a are consistent to within 
a — —2.30(4) = —0.383(6)^ except for a somewhat larger 
error at Jcnd — 0.4 and 0.6. (See Table U and Fig. 
[3l) Thus, we conclude that the value of a is indepen- 
dent of Jcnd- The estimated value of a is comparable 
to a = -0.32^ in Ref. In Fig. O the dotted line 

represents Jcnd = 0.50865 determined in Ref. f^]. 



IV. APPLICATION TO S = 2 HEISENBERG 
CHAIN AND CONCLUSIONS 

We have shown that the energy spectrum modified by 
the tuning parameter Jend can be fitted using an effec- 
tive massive relativistic dispersion with a boundary scat- 
tering length ab(Jend) modified for lattice models. The 
intermagnon scattering length a is constant irrespective 
of Jend, as wcU as other bulk quantities including the 
Haldane gap, the magnon velocity, and the correlation 
length. In contrast, ab(Jcnd) drastically changes around 
Jcnd ~ 0.5, representing a crossover point for the physics 
at the boundary. 
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TABLE I. Results of numerical simulations for a single itinerating magnon, showing magnon energy A, magnon velocity v, 
boundary scattering length a^, intermagnon scattering length a, and correlation length ^ — v/A, Ob/^ and a/^. 



J end 


A 


V 




a 











0.4104792487(1) 


2.466838(1) 


-0.3748(1) 


-2.30(2) 


6.009654(1) 


-0.06237(2) 


-0.383(4) 


0.1 


0.4104792486(1) 


2.466844(2) 


-0.0836(3) 


-2.301(2) 


6.009669(4) 


-0.01391(5) 


-0.3830(4) 


0.2 


0.4104792487(1) 


2.46684(1) 


0.540(2) 


-2.303(5) 


6.00966(3) 


0.0898(3) 


-0.3833(9) 


0.3 


0.4104792486(4) 


2.46684(4) 


2.081(8) 


-2.30(4) 


6.0096(1) 


0.346(1) 


-0.384(6) 


0.4 


0.410479248(2) 


2.4668(3) 


7.33(5) 


-2.3(2) 


6.0098(7) 


1.220(9) 


-0.38(5) 


0.6 


0.410479248(2) 


2.4668(2) 


-16.93(2) 


-2.3(2) 


6.0096(5) 


-2.821(6) 


-0.38(3) 


0.7 


0.4104792483(2) 


2.46685(3) 


-9.586(5) 


-2.30(3) 


6.00968(7) 


-1.5951(9) 


-0.382(4) 


0.8 


0.4104792483(2) 


2.46685(2) 


-7.317(3) 


-2.30(2) 


6.00968(4) 


-1.2176(5) 


-0.383(3) 


0.9 


0.4104792483(2) 


2.46685(2) 


-6.233(3) 


-2.30(2) 


6.00968(4) 


-1.0372(5) 


-0.383(3) 


1.0 


0.4104792485(1) 


2.46684(1) 


-5.605(3) 


-2.30(2) 


6.00967(4) 


-0.9328(5) 


-0.383(3) 







1 

o 


1 




1 


1 

Ob 
a 


o 

X 




(r 


o 


o 












; 


X 


X X 


X 


X X 


X 


X 














o 


o 












o 












1 


1 


o 
1 













0.2 


0.4 


0.6 


0.8 







end 



FIG. 3. (Color online) Boundary scattering length and inter- 
magnon scattering length as a function of Jcnd- The dotted 
line represents Jend = 0.50865. 



Analysis of the boundary scattering length and inter- 
magnon scattering length was also carried out for a. S = 2 
AF Heisenberg chain, where Si and in the Hamilto- 
nian in eq. ((IJ represent the S = 2 and 5 = 1 operators, 
respectively. In addition, we choose Jcnd = 1, so that 
the low- lying magnon states are IMs, and a similar for- 
mula for i?io is obtained to that shown in Eq. ([5]). Our 
data was taken using rris = 1024, which corresponds to 
rus" ~ 6000, and large systems up to L + 2 = 2048, 
The truncation error is smaller than 1 x 10^^^. Note 
that the numerical cost using NA-DMRG is about 200 
times less than that for the standard DMRG in this case. 
In contrast to a former report, our results suggests a 
large value of ab(Jo„d = 1) -33(1) = -0.67(2),^. Our 
calculations give the excitation gap A — 0.0891623(9), 
the spin velocity v = 4.42(1), and the correlation length 
^ = 49.6(1). In particular, the value of A obtained in the 
present study is consistent with the value of 0.08917(4) 



determined by quantum Monte Carlo simulations.'^'' The 
estimate of the Haldane gap has thus been improved by 
two more significant digits. This indicates the ability of 
the effective theory to correctly describe the low-energy 
physics, and the usefulness of the proposed numerical ap- 
proach is studying such problems 

It would be of interest to apply the approach used in 
this work to finite size scaling with different boundary 
tuning methods such as hyperbolic deformation4^i^i^. 
In such a situation, the excited quasi particle is weakly 
confined near the center of the system under the defor- 
mation. In Ref. [^^], we showed that it is necessary to 
introduce an additional parameter d and replace L + 1 
by L + d in order to reduce higher-order corrections. This 
replacement is introduced in the effective model shown in 
Eq. ([4]) by considering the effective boundary scattering. 
The boundary scattering length has an important and 
universal infiuence on excitation energy scaling as long 
as there are chain ends. 

In this work, a relation ^ = v/A is used to estimate 
the correlation length in each spin-S" chain. If we use 
an assumption for a relation between the low-energy dis- 
persion curve and the ground state correlation length, 
namely sinh^~^ = A/v in this case^^ the correlation 
lengths are evaluated as 6.03720(9) in S' = 1 and 49.6(1) 
in S — 2. In the case of S" = 1, we have a meaning- 
ful different value from the former estimation. On the 
contrary, the difference is not confirmed in the case of 
S = 2. To find correct relation between the low-energy 
dispersion and the correlation length in each spin-S" AF 
Heisenberg chain is a future issue. 

For a final development of the low-lying effective field 
theory to describe the low-lying magnon dispersions, dis- 
cussions for rigorous results of wave functions and energy 
dispersions for low- lying states are important. The effec- 
tive dispersion relation of \/ A^ -I- v"^ sin^ fceff is known to 
appear in the Haldane phased and also in the massive 
phase of the 5" = 1/2 XXZ model.^^ The Bethe-ansatz 
solutions for OBC suggest that an analogous crossover 
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from an IM with real fccg to a BM with a damping na- 
ture can be found as a continuous change from a real to 
an imaginary rapidity.— 
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